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The observed galactic magnetic fields may have a primordial origin. I briefly 
review the observations, their interpretation in terms of the dynamo theory, and 
the current limits on cosmological magnetic fields. Several possible mechanisms 
for generating a primordial magnetic field are then discussed. Turbulence and the 
evolution of the microscopic fields to macroscopic fields is described in terms of a 
shell model, which provides an approximation to the full magnetohydrodynamics 
and indicates the existence of an inverse cascade of magnetic energy. Cosmological 
seed fields roughly of the order of 10 -20 G at the scale of protogalaxy, as required 
by the dynamo explanation of galactic magnetic fields, seem rather plausible. 



1 Introduction 

Apart from the baryon number and the spectrum of energy density fluctua- 
tions, the physical processes that took place in the very early universe do not 
have many consequences that could still be directly detectable today. Most ob- 
servables have been washed away by the thermal bath of the pre-recombination 
era. One possibility, which has recently received increased attention, is offered 
by the large-scale magnetic fields^abserved in a number of galaxies, in galaxy 
halos, and in clusters of galaxies El. The astrophysical mechanism responsible 
for the origin of the galactic magnetic fields is not understood. Instead one 
postulates a small seed field, which can then be either enhanced by the com- 
pression of the protogalaxy, or exponentially amplified by the turbulent fluid 
motion as in the dynamo theorya. The exciting possibility is that the seed field 
could be truly primordial, in which case cosmic magnetic fields could provide 
direct information about the very early universe. 

The problem then is twofold. First, one has to find a mechanism in the 
early universe which is able to produce a magnetic field large enough to act as 
the seed field. There are various proposals, a number of which are based on 
the early cosmological phase transitions, which are discussed in Sect 3. The 
second problem is to explain how the initial field, which is expected to be 
random as it is created by microphysics and having correlation lengths typical 
to microphysics, can grow up to be coherent enought at large length scales. 
This is a problem in magnetohydrodynamics which is discussed in Sect 4. 
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2 Observations and limits 



2. 1 Observing cosmic magnetic fields 

Magnetic fields at the level of few fj,G have been detected in galaxies, in galactic 
halos, and in clusters of, galaxies. They can be observed indirectly at optical 
and radio wavelengths u. The Zeeman splitting of the spectral lines would 
provide a direct measure of the strength of the magnetic field, but the shifts 
are very small and this method is applicable mainly to our own galaxy. 

Electrons moving in a magnetic field emit synchrotron radiation, and both 
its intensity and polarization can be used to extract information about the 
magnetic field. However, one needs first to fix the relative magnitudes of the 
electron and magnetic field densities. Usually equipartition of magnetic and 
plasma energies is assumed. 

Information about distant magnetic fields in e.g. clusters in galaxies has 
been obtained by studying the Faraday rotation of polarized light. The method 
is based on the fact that the plane of polarization of linearly polarized elec- 
tromagnetic wave rotates as it passes through plasma supporting a magnetic 
field. The rotation angle Ax depends on the strength and extension of the 
magnetic field, the density of plasma, and on the wavelength A of radiation. 
The Faraday rotation measure (RM) is defined as 



where B\ is the strength of the magnetic field along the line of sight. This 
method requires some independent information about the electron density and 
the field reversal scale. 

2.2 a — lu dynamo 

The currently favoured explanation for the origin of the large scale galactic 
magnetic fields is the a — oj dynamo 13, which through turbulence and differ- 
ential rotation amplifies a small frozen-in seed field Bo to the observed /iG 
field. An initially toroidal seed field, which is carried along on the disc of a 
rotating (spiral) galaxy, is locally distorted into a loop by the up- or downward 
stochastic drift of the gas. As the gas moves away from the plane of the disc, 
the pressure decreases and the gas expands; at the same time it is subject to 
a Coriolis force which will rotate it. The magnetic field lines, glued to the gas, 
will follow and thus a poloidal component perpendicular to Bo is generated. 
The small poloidal loops will reconnect and coalesce to produce a large scale 
field. This is the so-called a-effect. Because the disc does not rotate like a 
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rigid body, the field lines will be wrapped and the poloidal field will induce a 
toroidal field; this is the w-effect. Thus the dynamo mechanism does not only 
produce a large scale field but can also, to some extent, predict the shape of the 
field. In practise one does not attempt to describe the small scale turbulence 
directly but rather resorts to a mean-field dynamo, where the average induced 
current is given by j,- n( j = eraBo, where a is conductivity and a is a parameter 
related to the average drift velocity. 

The dynamo saturates when the growth enters the non-linear regime. A 
typical growth time is of the order of 10 9 years, with the rotation period of 
the galaxy setting a lower limit on the growth time. Recently it has been 
argued El that the saturation might actually be too fast for a large-scale field 
to form. The strength of the required initial seed field is rather uncertain, but 
as a rule of thumb one could use a value like 10 -20 G on a comoving scale of 
a protogalaxy (100 kpc). 

Another possibility is that the galactic field results directly from a pri- 
mordial field, which gets compressed when the protogalactic cloud collapses. 
The primordial field strengths needed are however quite large. In any case, 
it appears as if a primordial field is required to explain the observed galactic 
magnetic fields (although it is conceivable that a purely astrophysical solution 
could exist as well). 

2.3 Limits on cosmologial fields 

Observing cosmological, intergalactic magnetic fields would be of great impor- 
tance and would strengthen the case for their truly primordial origin. It has 
been claimed by PlagaQ that the arrival times of extragalactic 7-rays could be 
used to detect fields as weak as 10 -24 G. The idea is that cosmic rays, orig- 
inating from far-away objects such as QSOs or gamma-ray bursters, scatter 
off the background cosmic magnetic field. This gives rise to pair production 
and a delayed 7-ray which could then be observed, and the ratio of prompt to 
delayed 7's provides a measure of the strength of the intergalactic magnetic 
field. The energy spectrum of 7-rays is likewise affected by the intergalactic 
magnetic fieldtl. Ultra-high cosmic ray protons would also be deflected, by the 
intergalactic magnetic fields so that their arrival times could be usedH to set 
bounds as low as 10~ 12 G. 

A direct limit on cosmological magnetic field has been obtained by consid- 
ering the .rotation measure of a sample of galaxies and quasars, which yields 
the limit B < 10 -9 G (fiichioo/0.01)- 1 , where ft IG is the fraction of the 
ionized gas density of the critical density in the intergalactic medium, and h 
is the Hubble constant. It has also been suggested that the power spectrum of 
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the cosmological field could he determined by studying the correlations in the 
RM for extragalactic sources a. 

If a primordial magnetic field is present at the time of recombination, it 
will give rise to anisotropic pressures^. These would distort the microwave back- 
ground. Barrow, Ferreira and Silk El have recently considered the effect of a 
magnetic field on the evolution of shear anisotropy in a general anistropic flat 
universe. They then compared the result with the 4-year COBE data set, and 
assuming that the whole observed anisotropy is due to magnetic stresses, con- 
cluded that if the field can be taken homogeneous, B < 3.4xlD -9 G (f^ohlo) 1 / 2 . 
This is a more stringent limit than what can be obtained Ej from primordial 
nucleosynthesis considerations. These take into account the contribution of the 
magnetic field to the Hubble expansions rate and to the n <-> p reaction rates, 
plus the fact that the electron wave functions get modified (the phase space of 
the lowest Landau level is enhanced). The limit thus obtained is B < 10 11 G 
for a constant B at T = 10 9 K, and if the field is inhomogeneous, the limit is 
less stringent by an order of magnitude. 

The anisotropy limit also indicates that the distortions of the Doppler 
peaks of the microwave background til due to a homogeneous magnetic field 
are unobservable. The Faraday rotation in the polarization of the microwave 
background could still be observablelia. 

3 Generating primordial magnetic fields 

There are a number of proposals for the origin of magnetic fields in the early 
universe, many of them involving the early phase transitions. Fluctuations in 
the electromagnetic field in a relativistic plasma .axe by themselves sufficient 
for generating a small scale random magnetic fieldtLil To obtain a seed field for 
galactic magnetic fields, one however needs fields with much larger coherence 
lengths. 

3.1 Spontaneous charge separation 

Magnetic fields may arise in a electrically neutral plasma if local charge separa- 
tion happens to take place, thus creating a local current. Itiias_been proposed 
that this could occur during a first order QCDEJ or EWEJ'EEl phase transi- 
tion, which proceed by nucleating bubbles of the new phase in the background 
of the old phase. There one finds net baryon number gradients at the phase 
boundaries, providing the basis for charge separation, and the seed fields arise 
through instabilities in the fluid flow. For that one has to require that the 
growing modes are not damped. Turbulent flow near the walls of the bubbles 
is then expected to amplify and freeze the transient seed field. The various 
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hydrodynamical features have been studied in linear perturbation theory by 
Sigl, Olinto and JedamzikliS, who argued that on a 10 Mpc comoving scale, 
field strengths of the order of 10~ 29 G for EW and 1CT 20 G for QCD could be 
obtained. These might furthex.be enhanced by several orders of magnitude by 
hydromagnetic turbulenceEZrtj, which will be discussed in Section 4.3. 



3.2 Bubble collisions 

In a first order phase transition the phases of the complex order parameter 
<J> = pe l6 /^/2 of the nucleated bubbles are uncorrelated. When the bubbles 
collide, there arises a phase gradient which acts as a source for gauge fields. 
The phase itself is noi^a gauge invariant concept, but it has been pointed out by 
Kibble and Vilenkinlij that a gauge invariant phase difference can be defined 
in terms of an integral over the gradient D^Q. 

Magnetic field generation in the collisianxif phase transition bubbles has 
been considered in the abelian Higgs modelEflo. One assumes that inside the 
bubble the radial mode p settles rapidly to its equilibrium value r\ and can thus 
be treated as a constant. The dynamical variables are thus O and the gauge 
field A^. The starting point is the U(l)-symmetric lagrangian 

C = --F^ v + D^(D^ + V(\$\), (2) 

where the potential V is assumed to have minima at p — and p = r/. The sim- 
plest case is that two spherical bubbles nucleate, one at (x, y, z, t) — (0, 0, 0) 
and the other at {x,y,z,t) ~ (0, 0, — zq, to)j an d keep expanding with the ve- 
locity v even after colliding. Because of the symmetry of the problem, the 
solutions to the equations of motion are functions of z and r = \Jt 2 — x 2 — y 2 
only, and the U(l) gauge field = x M /(r, z). With the appropriate initial 
conditions, the solutions arelij 
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where to = \Jk 2 + e 2 r) 2 , R is the radius of the bubbles at the collision and Z\ 
the point of first collision on the z-axis. Here the velocity of the bubbles has 
been taken to be v = 1. The generated magnetic field is rapidly oscillating 
and orthogonal to the z-axis and in this case confined inside the intersection 
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Figure 1: The time evolution of B along the radius r = 1 for t=10, 20, 30 and 40 after the 
initial collision. The radius of the bubbles at the collision has been chosen here R = 10, and 
the collision point on the z-axis is z\ = 50. (The units are such that er\ = 1). 



region, as can be seen from Fig. 1. It has a ring-like shape in the (x,y)-planc. 
It can be shown that subsequent collisions of the bubbles, which now may have 
a magnetic held inside the bubbles, nevertheless lead to a qualitatively similar 
outcomeEI 

There are however two additional important ingredients which need to be 
taken into account: the high but finite conductivity E3t3 of the primordial 
plasma and the fact that in the electroweak phase transition the bubbles will 
in fact intersect with non-relativistic velocitiescJ. Finite conductivity gives rise 
to diffusion, the consequence of which is to smooth out the rapid oscillations 
of B, whereas low v will permit the magnetic flux to escape the intersection 
region and penetrate the colliding bubbles, where its evolution will be governed 
by usual magnetohydrodynamics. The resulting behaviour is demonstrated in 
Fig 2. 

The strength of the generated magnetic field depends on the bubble wall 
velocity in an essential wayEj. In the electroweak case the initial growth of the 
bubble wall is by subsonic deflagration, with velocities of the order of 0.05c, 
depending on the assumed friction strengthen. The wall is preceded by a shock 
front, which may collide with the other bubbles. This results in reheating, and 
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Figure 2: B along r = 1 with R = 10 at t = 30 after collision with v = 0.8 and 0.1. Here 
conductivity has been taken to be a = 7T The point of initial collision on the z-axis is 
2i = 50, and the outer edge of the intersection region is at ~ 50 ± vt. The outer edge of the 
magnetic field is at ~ 50 ± t (and units are er) = 1). 



oscillations of the bubble radii, but eventually a phase equilibrium is attained. 
The ensuing bubble growth is very slow and takes place because of the ex- 
pansion of the universe. Because the universe has been reheated back to T c , 
no new bubbles are formed during the slow growth phase. Assuming that the 
abelian Higga model results are applicable to the electroweak case, as seems to 
be the caselPJ, one may estimate thatEJ 



B ~ 2.0 x 10 2 V7 2 + 2-yR/R G, (4) 

where it has been assumed 0o — 1, T c = er) = 100 GeV; the average dis- 
tance between the nucleation centers r ave = 9.5 x 10 -8 £# and the velocity 
V = 1.2 x 10~ 4 were taken as reference values. Folding in the spectrum of 
separation of the adjacent shocked spherical bubblesc3, averaging over all pos- 
sible inclinations of the ring-like magnetic field, and taking into account the 
enhancement of magnetic energy due to an inverse cascade (see next Section), 
one arrives at the estimate B—s ~ 10~ 21 G for the cosmological magnetic field 
at the scale of 10 Mpc today E3. 



3.3 Fluctuating Biggs gradients 

It has been pointed out by VachaspatiB that fluctuating Higgs field gradients 
will induce a magnetic field. Such local fluctuations are naturally present at 
the EW phase transition, since the embedding of the electromagnetic field in 
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SU(2) ® U(1)y involves these gradients: 

Ft™ = -WjVj-ViV}) , 

v > - w/¥^ ■ 

where <fi is the Higgs field. (An_argument against horizon-size fluctuations 
has been presented by Davidson EH) . At the electroweak phase transition the 
correlation length in the broken phase is ~ 1/mw (assuming that the Higgs 
mass is comparable to %). The field strength Ff™ is thus constant over a 
distance ~ 1/mw, but it varies in a random way over larger distances in order 
to respect causality. The vector Vi is also random, of course. Its variation is 
due to the fact that the Higgs field 4> makes a random walk on the vacuum 
manifold of <j>. The problem then is to estimate the field strength FfJ 71 over a 
length scale ~ N/mw- If N = 1, then it follows on dimensional grounds that 
the root-mean-square F^ 71 ~ m w ~ 10 24 G. For N large, one should use a 

statistical argument. Vachaspati argued that the gradients are of order 1/y/N, 
since <f> makes a random walk on the vacuum manifold with A</> ~ y/N, and 
since Ax ~ N. Thus Vi is, in a root mean square sense, of the order 1/y/N, 
and hence Ff" 1 is of order 1/N. Taking further into account that the flux in 
a co- moving circular contour is constant, the field must decrease like l/a(t) 2 , 
where a(t) is the scale factor. However, there also exists a different statistical 
scenario, based on a line average, where the gradient vectors are taken to be the 
basic stochastic variablescS. The field strength can be interpreted statistically 
in such a way that the mean magnetic field satisfies 

(F-)t =0 , . (6) 

One observes that the scaling behavior is weaker by a factor y/N. This means 
that for a scale of 100 kpc 

yJ(F?J todav ~ 10- 18 G , (7) 

which is very close to the value desired for the dynamo effect. Because the 
Faraday rotation involves the average along the line of sight, this averaging 
method could indeed be appropriate for those observations. 

3-4 Vacuum condensates 

Another, more exotic possibility for generating primordial magnetic fields is 
based on the observation that, due to quantum fluctuations, the Yang-Mills 
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vacuum is unstable in a large enough background magnetic field at zero 
temperature. There_are indications from lattice calculations that this is a non- 
perturbative result EJ. In a pure SU(N) theory at the one-loop level the zero 
temperatum effective energy for a constant background non-abelian magnetic 
field reads 



V(B) = \b 2 



117V 
96tt 2 



9 2 B 2 



li 2 2 



with a minimum at 
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A8tt 2 
llNg 2 



(8) 



(9) 



and Knin = V(B m [ n ) = — 0.029(g£? m i n ) 2 • Thus the ground state has a non-zero 
non-abelian magnetic field, the magnitude of which is exponentially suppressed 
relative to the renormalization scale, or the typical momentum scale of the sys- 
tem. In the early universe, however, where possibly a grand unified symmetry 
is valid, the exponential suppression may be less severe. It is also attenuated 
by the running of the coupling constant. 

In the early universe the effective energy picks up thermal corrections from 
the fermion, gauge boson, and Higgs boson loops. The detailed form of the 
thermal correction depends on the actual model, but one may take the cue 
from the SU(2) one-loop calculation, where they are obtained by summing the 
Boltzmann factors exp(— j3E n ) for the oscillator modes 



E 2 



p 2 +2gB(n+-] 



2gBS 3 + m 2 (T), 



(10) 



where S3 = ±1/2 (±1) for fermions (vectors bosons). Eq. ( |l0| ) includes the 
thermally induced mass m(T) ~ gT, corresponding to a ring summation of the 
relevant diagrams. Numerically, the effect of the thermal mass turns out to be 
very important. At high temperature, the leading behaviour is gip£en mainly 
by the bosonic contributions, and thus one arrives at the estimate Ej 
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(11) 



where Kf{x) = gBl 2 /(AxT 2 ) + m 2 a x/(gB) and numerically SVf ~ 0.02 (gB) 2 
which serves only to shift the value of B at the vacuum slightly. Thus B ^ 
is the state of the lowest energy even at high temperature. 

A local fluctuation will then trigger the creation of the a new vacuum 
with non-zero non-abelian magnetic field inside a given particle horizon at 
scales jtx ~ T. The Maxwell magnetic field is then just a projection of the 
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non-abelian field, the order of magnitude of which is given by Eq. (^). One 
can then estimate that EJ 

Bm^exp^X^Sx^^) 2 G, (12 ) 

where the reference number is for susy SU(5). 

Electrxiweak magnetic condensates have also recently been considered by 
Cornwall E3, who has suggested that they would give rise to magnetic fields 
with a net helicity via_the generation of electroweak Chern-Simons number. 
Joyce and ShaposnikovEj have argued that a right-handed electron asymmetry 
generated at the GUT scale, could give induce a hypercharge magnetic field 
via a Chern-Simons term. Both suggestions deserve further study. 



3. 5 Inflation and magnetic fields 

The basic problem with inflation with regards to magnetic field generation is 
that the early universe was a good conductor so that, ignoring turbulence, 
the magnetic flux ~ Ba 2 tends to be conserved, where a is the cosmic scale 
factor. To avoid this, one needs to break the jxmformal invariance somehow, 
as was first suggested by Turner and WidrowEJ, who considered couplings to 
the curvature R such as RF 2 and RA 2 , as well as photop^axion couplings. 
Dilaton coupling of the form e®F 2 has also been considered^.], and interesting 
field strengths can be obtained at the expense of tuning the coupling strength. 
If a phase transition takes place during the inflationary period, a sufficiently 
large magnetic field can be created, provided however that the phase transition 
takes place during the final 5 e- foldings E3. 



4 Prom microscopic to macroscopic 

Even assuming that a primordial magnetic field is created at some very early 
epoch, a number of issues remain to be worked out before one can say anything 
definite about the role of primordial fields for the origin of galactic magnetic 
fields. At the earliest times magnetic fields are generated by particle physics 
processes with length scales typical to particle physics. Jit has been shown 
that such fields are stable against thermal fluctuations E3.) The remaining 
question is whether it is at all possible for the small scale fluctuations to 
grow to large scales, and what exactly is the scaling behaviour of B ms or the 
correlator (B(r + x)B(x)). To study these problems one needs to consider the 
detailed evolution of the magnetic field to account for such issues as to what 
happens when uncorrelated field regions come into contact with each other 
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during the course of the expansion of the universe. In general, turbulence is 
an essential feature of such phenomena. These questions can only be answered 
by considering magnetohydrodynamics (MHD) in an expanding universeEH. 



4-1 MHD in curved space 

Let us consider the early universe as consisting of ideal fluid with an equation 
of state of the form p — ^p, where p is pressure and p the energy density. 
Let us further assume that the fluid supports a (random) magnetic field. The 
energy-momentum tensor is then given by 



= (p + p)U^U v +pg^ 

+ ^(F^F» a -±g^F Xa F x °y (13) 

and is the four- velocity of the plasma, normalized as U^U^ = —1, and 
Fu V = duA v — d u A fjL is the electromagnetic field tensor. Note that, as long as 
diffusion can be neglected, the presence of the magnetic field does not change 
the equation of state. The magnetic energy is further assumed to be much 
smaller than the radiation energy, so that one can assume a flat, isotropic and 
homogeneous universe with a Robertson- Walker metric. Although the mag- 
netic field generates local bulk motion, this may still be consistent with isotropy 
and homogeneity at sufficiently large scales, in particular if the magnetic field 
is random, i.e. statistically homogeneous and isotropic on scales much larger 
than the intrinsic correlation scale of the field. 

For numerical treatment it is convenient to write the equations of mo- 
tion explicitly in 3+1 dimensions and use the conformal time_t = J dt/a as a 
variable. After a lengthy derivation one arrives at the formsEJ 

— = —(V • v)S — (v • V)S — Vp + J x B. (14) 



and 



dt 

2 7 2 + l dlnp _ dS 2 /di 



4 7 2(2 7 2-l) dt (| / 5 7 ) 2 (2 7 2-l) 



v Vln(p 7 2 )- V-v+i^, (15) 



4 ~ ? ' 



where we have set p+p = -*p. The Maxwell equations can be written explicitly 



3 

as 



— =-VxE, V-B = 0, (16) 
at 
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and 

J = VxB-- T , V-E = p e (17) 
at 

where p e is the charge density and p e = a 3 p e . Further, 

E = -vxB, (18) 

which is valid in the limit of high conductivity. 

The evolution of a random magnetic field configuration in 2d is shown in 
Fig. 3 for a lower and higher resolution. As time goes on, the coalescence 
of magnetic structures is seen to lead to the gradual formation of larger and 
larger scales. Such a behaviour is encouraging, indicating that the field is not 
really comovingly frozen. On the other hand, 2d MHD is very different from 
the 3d MHD, and moreover, the Reynolds number in the-.simulation (about 
10) is wildly unrealistic. However, Dimopoulos and DavisEa have also pointed 
out that when two initially uncorrelated domains come into contact, the field 
at the interface should untangle with the plasma bulk velocity v to avoid the 
creation of domain walls. They propose that the correlation length £ evolves 
according to 

+ (19) 
where H is the Hubble parameter and the velocity v depends dynamically 



on B and should, in principle, be determined from MHD. Nevertheless, (19) 
again points towards the possibility of the magnetic field not necessarily being 
comovingly frozen. 



4-2 Shell models 

In ordinary hydrodynamics many properties of turbulence, in particular those 
related to energy transfer and to the spectral properties have been studied 
successfully using a simple cascade model. This is true not only qualitatively, 
but also quantitatively, which is the reason why the cascade model is now much 
used in studies of nonlinear physics lj. 

The basic idea is that the interactions due to the nonlinear terms in the 
MHD equations are local in wavenumber space, and in fc-space the quadratic 
nonlinear terms become a convolution. Interactions in fc-space involving tri- 
angles with similar side lengths have the largest contribution. This has led to 
the shell model which is formulated in the space of the modulus of the wave 
numbers. This space is approximated by N shells, where each shell consists 
of wave numbers with 2" < k < 2™ +1 (in the appropriate units). The Fourier 
transform of the velocity over a length scale k~ x (k n = 2") is given by the 
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Figure 3: Left column: magnetic field lines at different times at low resolution (64 X 64 
meshpoints). Right column: magnetic field lines at different times at higher resolution 
(128 X 128 meshpoints). 
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complex quantity v n , and B n denotes a similar quantity for the i3-ficld. Fur- 
thermore, the convolution is approximated by a sum over the nearest and the 
next nearest neighbours, 

2 

N n (v,B) = C l0 v n+t B n+J . (20) 

i,J=-2 

Here v and B have lost their vectorial character, which reflects the fact that this 
model is not supposed to be an approximation of the original equations, but 
should be considered as a toy model that has similar conservation properties 
as the original equations. 

Velocity and magnetic fields are thus represented by scalars at the discrete 
wave numbers k n = 2" (n = 1, N), i.e. k n increases exponentially. There- 
fore such a model can cover a large range of length scales (typically up to 
ten orders of magnitude) . The important conserved quantity is i?tot-R 4 , whexe 
Etot = J T 00 d 3 x is the total energy. The resulting equations of motion readliU 

lp ^ = Nn (v,b), (21) 

^ = M n (v,b), (22) 
at 

where 

2N n (v, b) = ik n (A + C)« +1 < +2 - b* n+1 b* n+2 ) 

+ik n (B - - b* n _ x b* n+x ) (23) 

-ik n {\B + |A)«_ 2 <-i - ^-2^-1). 

M n (v, b) = ik n (A - C)(v* n+1 b* n+2 - b* n+1 v* n+2 ) 

+ik n (B + iCXiCAVi " K-i< + i) (24) 
-ik n {\B ~ |^)«_ 2 fo*_ 1 - K-2<-l)> 

with A, B, and C being free parameters. It is straightforward to verify that 
2 E < N n + E K M n = 0, using that k n = 2". 

4-3 Inverse cascade 

The numerical study of the cascade model requires of course that the param- 
eters A, B, C are fixed so that the model has the same conservation laws as 
the full-fledged MHD. The model can then be solved numerically Ell, and the 
results are shown in Fig. 4, where the transfer of magnetic energy to larger and 
larger length scales is clearly seen. This process, the inverse cascade, is due 
to the nonlinear terms giving rise to mode interactions. The initial magnetic 
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energy spectrum was chosen to be given by Em(t — 0) ~ k with the total 
magnetic energy equal to pq. The number of shells was N — 30 so that length 
scales differing by ten orders of magnitude were covered. It was found that 
the integral scale, which measures where most magnetic energy is concentrated 
and which is given by 

l = J {2ir/k)E M (k)dk I J E M (k)dk, (25) 

where EM(k) is the magnetic energy spectrum, increases with the Hubble time 
approximately like t^ 25 . 

However, around the time of recombination the photon mean free path A 7 
became very large and photon diffusion became very efficieniin smoothing out 
virtually all inhomogeneities of the photon-baryon plasma EJ. This process is 
often referred to as Silk damping, which corresponds to a kinematic viscosity 
v ~ A 7 (in natural units). Silk damping may thus destrxty the magnetic field, 
as has been noted by Jedamzik, Katalinic and OlintoEj. Therefore one has 
to follow numerically the evolution of the magnetic and kinetic energy spectra 
in the presence of kinematic viscosity. The results E3 are presented in Fig. 
5 and the main point can be summarized as follows: in the cascade models 
magnetic energy is transferred to large length scales even in the presence of 
large viscosity. Here the initial magnetic spectrum was chosen to be flat in 
accordance with the large time behaviour suggested by Fig. 4. 

Foii-a sufficiently large viscosity, the inverse cascade stops. One may es- 
timate tLS that this typically takes place close to recombination. The results 
suggest that in the real MHD, inverse cascade is operative and is essentially 
not affected by Silk damping, except very late and perhaps for very weak fields. 
Thus we may conclude that it is unlikely that an equipartition exists in the very 
early universe. A similar conclusion can be drawn in a different, continuous 
model where the,inverse cascade can be found analytically in an appropriate 
scaling regimetL 



5 Conclusions 

Explaining the galactic magnetic fields in terms of microphysical processes 
that took place when the universe was only ten billionth of a second old is 
a daunting task, which is not made easier by the complicated evolution of 
the magnetic field as it is twisted and tangled by the flow of plasma. It is 
nevertheless encouraging that mechanisms for generating primordial magnetic 
fields of suitable size exist, and in particular those based on the early cosmo- 
logical phase transitions discussed in Sect. 3 look promising. At the same 
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Figure 4: Spectra of the magnetic energy at different times. The straight dotted-dashed 
line gives the initial condition (to = 1), the solid line gives the final time (t = 3 X 10 4 ), and 
the dotted curves are for intermediate times (in uniform intervals of Alog(f — to) = 0.6). 
A = 1, B = -1/2, and C = 0. 
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Figure 5: Magnetic and plasma kinetic energy spectrum as a function of the wave number k 
in the cascade model for small (y = 10~ 2 ) plasma viscosity (left) and large (v = 10 2 ) plasma 
viscosity (right). The highest time is t = 10 8 , corresponding to a Hubble time 10 16 . 
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time the fact that there are so many possibilities tend to underline our igno- 
rance of the details of the subsequent evolution of the magnetic field. The 
step from microphysics to macroscopic fields is a difficult one because of the 
very large magnetic Reynolds number of the early universe. However, different 
considerations, both analytic approximations, 2d simulations, as well as the 
full-fledged shell model computations which can account for turbulence, seem 
to point to the existence of an inverse cascade of magnetic energy. Moreover, 
as discussed in Sect. 4.3, the inverse cascade is obtained also in the presence 
of a large plasma viscosity. Therefore the primordial origin of the galactig 
magnetic fields is quite possible. 

Much theoretical work remains to be done, though. At the same time it is 
very important that progress is made on the observational front. In particular, 
measuring or setting a stringent limit on the intergalactic field, which could be 
possible in the near future as indicated in Sect. 2.3, would provide the testing 
ground for all theoretical scenarios. 
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